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Abstract 


For the prediction of fan noise, the two- 
dimensional method has been used in many cases. 
However, the accuracy of the two-dimensional method 
needs to be assessed against three-dimensional 
effects, Therefore, two-dimensional, quasi three- 
dimensional and three dimensional theories for the 
prediction of pure tone fan noise due to the inter- 
action of inflow distortion with a subsonic annular 
blade row have been studied with the aid of an un- 
steady three-dimensional lifting surface theory 
developed by Namba. 1 ) 2 The present study also 
examines the effects of compact and noncompact 
source distributions on pure tone fan noise in an 
annular cascade. Numerical results show that the 
strip theory and quasL three-dimensional theory are 
reasonably adequate for fan noise prediction. The 
quasi thrco-dimonslonol method is more accurate for 
acoustic power and modal structure prediction with 
on acoustic power estimation error of -about ±2 dB. 
This accuracy seems to bo affected not by the sub- 
rcsonancc or super-resonance status of the acoustic 
modes, but by the reduced frequency of the inflow- 
distortion and by the inflow distortion radial dis- 
tribution. Also, the compact source prediction is 
different from the noncompnct source prediction by 
as much as 15 dB in upstream radiation cases and by 
oo much as 6 dB in downstream cases, depending bn 
reduced frequency and intorbladc phase. 


Introduction 


A major source of fan noise is unsteady flows 
Interacting with blades and vanes, These unsteady 
flows ore inlet distortions- 1 " 2 such os .ltmoBphertc 
turbulence, ® cross-wind effects, 5 static test in- 
stallation wakes, fan In Let boundary layers mid the 
ground vortices that are apparent during many 
static test conditions, all of which can produce a 
significant levels of noise with a fan nr a turbo- 
fan engine. Because these inf Low conditions can 
vary during the course of a test program, and can 
vary significantly between static testing and for- 
ward flight of an engine, 10 the evaluation of noise 
reduction features using static testing and die 
extrapolation to flight of static noise levels is 
extremely complicated and may prove erroneous. 


For flyover noise prediction, source noise 
error is only n part of the problem since atmo- 
spheric propogation effects and installation effects 
also appear to contribute significantly. However, 
the study of fan noise sources Is an important part 
of this problem, and knowledge of the importance 
and strengths of possible sources of unsteady flows 
described above is particularly important in sel- 
ecting a facility for simulating flight behavior 
statically. The degree of Importonce of these un- 
steady flows will depend on the particular fan 
stage and unsteady flow structure. It is therefore 
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lraportont; to correlate the results olf various fan 
noise calculations with the existing experimental 
model data end thereby define the Interplay of 
Scales and intensities in determining noise levels 
and to suggest the most fruitful inflow control 
strategics. These calculations can also provide 
information about fan noise modal structure for 
acoustic suppressor design. 11 

Many theoretical studios af fan noise have 
boon based on two-dimensional linear cascade 
madels 2t “®i L2-14 that have made clear Che fundament- 
al features of the sound generation process and 
hove shown the effect of the various factors that 
influence sound source noise. However, it has not 
been made clear over what range of conditions these 
current two-dimensionul theories are able to pre- 
dict the unsteady blade forces and the sound field 
for the case of a three-dimensional annular blade 
row and for complicated inflow-distortions. The 
two-dimensional problem Is easier to handle and to 
obtain numerical results for; however, comparison 
of two and three-dimensional results is needed over 
a wide range af parameters and Inflow conditions, 

The available theoretical results on unsteady 
blade forces for a single airfoil in an oblique 
gusc 1 ^ allow that the unsteady lift force amplitude 
decrcuscs as the spanwisc wave number of the gust 
Increases. We can infer then, that three- 
dimensional effects will play an important role on 
fluctuating blade forces in the caso of interaction 
witli on inflow distortion of large radial extent. 

As Tyler and Sofrin 1,) pointed out, the interaction 
field generates many modes with various circumfer- 
ential and radial mode numbers, and whether these 
modes propagates or decay depends strongly on the 
three-dimensional duct cutoff phenomena. Therefore, 
Che role of three-dimensional effects in the pre- 
diction of unsteady blade forces and acoustic power 
and acoustic modal structure lias to be clarified. 

A theoretical analysis of the unsteady force 
and the total acoustic power in a three-dimensional 
annular cascade was carried out by Namba. 1 Several 
ot the three-dimensional effects oil unsteady forces 
and ncousLic power were obtained. However, In his 
numerical calculations, the reduced frequency is 
treated as a constant vaLue, without regard to in- 
flow distortion circumferential mode number and the 
rotational speed of rotor, and so correct numerical 
results for rotor-inflow distortion interaction 
noise are not presented, Also, in his calculation, 
the circumferential arid radial distributions of the 
inflow distortion were limited to the case of sim- 
ple harmonic distortion. 

Therefore, in this paper, a theoretical analy- 
sis of pure tone fan noise generated by Inflow 
distortion-rotor interaction was carried out with 
the aid of the three-dimensional unsteady lifting 
surface theory developed by Namba, 1 ) 2 accepting as 
input various inflow distortions, in order to pre- 
dict the forward 'and aft radiated pure cone energy 
and the modal energy distribution. Special empha- 
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sis Is placed upon the clarification of the accu- 
racy of available two-dimensional theory for the 
prediction of pure cone fan noise duo to the Inter- 
action of inflow distortion with a subsonic annular 
blade row. For this study, three-dimensional cal- 
culations, two-dimensional calculations (3trlp 
theory) and quasi three-dimensional calculations 
are carried out and compared. In the quasi three- 
dlmenstonol calculation, two-dimensional chorduisc 
dipole distributions at several radial positions 
are resolved into Fourior-Bessel dipole distribu- 
tions on an annular blade row, and then, with this 
dipole distribution, the pure tone fan energy and 
modal structure in a three-dimensional annular duct 
arc calculated. In addition, the effects of com- 
pact and noncompact sources in an annular blade row 
on pure tone fan noise arc also studied. 

The procedure and assumptions introduced In 
the present study are described as follows. The 
theoretical model consists of a 3lnglc three- 
dimensional annular cascade rotating f*t constant 
angular velocity in an annular rigid-wall duct of 
infinite axial extent. Thus, the duct end reflec- 
tion, the effect of upstream or downstream blade 
rows, and the effects of duct area variation are 
not considered. These effects except for the last 
one, could be Included in the solution procedure of 
the present study." The undisturbed flow Is uniform 
axially and the relative velocity over the whole 
blad'j span is limited to the subsonic range. In 
ordf/t to make the problem tractable, linearized 
theory Is adopted; that Is, not only fluctuating 
quantities due to cound but also those due to the 
inflow distortion converted with the basic flow, 
arc assumed to be small quantities of the first 
order in comparison to the mean pressure or the 
main stream velocity, respectively. 

Tile calculation procedure Is as follows. 

First, the periodic velocity fluctuations ut the 
rotor blades due to the inflow distortion are ob- 
tained, Second, In order to cancel the fluctuating 
component ol the velocity normal to the airfoLls, 
acoustic dipoles are assumed to be distributed over 
the surface of the rotor blade. Fluctuating pres- 
sure and velocity fields induced by the rotor 
blades are expressed through a kernel function in 
terms of tiie dipole distribution on a blade. The 
cpnnwlsc distribution of acoustic dipoles is given 
as a sum uf acoustic radial mode components and the 
kernel function Is ’resolved Into the corresponding 
modal compiSnente. Finally, the boundary condition 
that the fluctuating velocity normal to the aero- 
foil should vanish at each blade surface position 
Is introduced, Ily solving the resulting integral 
equation numerically, we obtained the dipole dis- 
tributions, which arc then used to calculate the 
.sound pressure In the far field. Thus, the pure 
tone acoustic power and modal structure upstream 
and downstream of the rotor are obtained. Tire in- 
’ f low distortion Is resolved into Its Fourier-Bessel 
components, and then the effects of each inflow 
.distortion component on the acoustic power Is super- 
imposed to obtain the total pure tone acoustic 
power. It is worth noting that the same acoustic 
mode can be generated by all the inflow distortion 
components that have the same circumferential mode 
number, oven though their radial mode numbers are 
different. 
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Analytical Model 

F luctuating Pressure Induced^ by Rotor Blades Row 

Hie theoretical model considered here consists 
of o single annular blade row with N blades ro- 
tating at a constant angular velocity In an 

annular rigid walled duct of infinite axial extent 
(Fig. 1). The fluid flow la composed of an uiw 
disturbed flow with a uniform axial velocity W a 
and small fluctuating flows due to the rotor blades 
and Inflow distortion. The flow Is inviscld, of 
uniform entropy and has no thermal conductivity. 

The fluctuation induced by the rotor is assumed to 
be isentroplc. The fluctuations of the fluf.d flow 
are small compared with the undisturbed flow. It 
is also assumed that the .fluid velocity relative to 
the blades Is subsonic along the whole spon and 
that the blades have no steady load, 

A cylindrical polar coordinate system Is used 
with axes (r,0,z) fixed to the rotor as shown In 
Fig. 1. The radial and axial coordinates r and 
z are normalized with respect tp the radius r£ 
at the blade tip. The dimensionless time £ is 
correspondingly scaled with respect to r$/W a . 

Then in accordance with the assumptions given 
above, the continuity, momentum and energy equa- 
tions are combined to give the following linearized 
equation 2 • *■' for the fluctuating pressure p 



'fills equutiou Is nothing other than the acoustic 
wave equation in the rotor-fixed coordinate and Its 
solution must Ballsty the boundary condition re- 
quiring that tile rad la L velocity vanish on the 
rigid walL surfaces at the hub and tip radii. In 
the present approximation, a rigid, thin nirfoll 
plnceJ in the. periodically fluctuating velocity 
field is directLy equivalent to a sheet of acoustic 
dipoles, or, in other words, an unsteady lifting 
surface. Therefore, a rotor blade row can be re- 
presented by N acoustic dipole sheets with the 
dipole axes normal to the surface. The N acoustic 
dipole sheets have a constant intcrblade phase dif- 
ference and are circumferentially equally spaced at 
the locations (p, <p + 2jir/N, ()(J - 0, 1, 2, . . ., 
N - l) , The fluctuating pressure due to the rotor 
blade row cun be expressed in the form 


p(r»0.*) = Yj ; X) , p • z ) ciux: . (3) 

qo -CO pa-M 


where 
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where C n j,(£,q,p) and F n ,i(;,q,p) arc the blade 
force coefficients defined by 


L P n>l«»‘ , * p) J 4» 


*4. 1 fr 


dp 


yp-tn JbQ 


l/pj 

( 10 ) 


X exp jin (0 -qiJ + lCmhp+UJKz • 0^/4 “ n „,ll z 'U] 


(7) 


let 


£P 0 ,p(p,t;)*e denotes the pressure jump across 
the blade surface caused by (q,p) mode components 
of the inflow distortion, where q and p are, 
respectively, a circumferential mode number and a 
radial mode number of the. Inflow distortion. 

C 0 denotes the axLal chord length of the blades, 
which is assumed^ constant along the span. 
ttC a *_ct (*l^r) r^/Wfl is the reduced frequency, 

M„ - W c /a 0 is the axial Mach number, PS ° 1 - Mg 
is the frond tl-Glauert parameter and further On, l 
Is defined by 


O' 


n, l 


fc.i - (» + vWJQfc < 8 > 


B n (kn,l‘r) and k n r nlc respectively the normal- 
lied radial eigen-function and its eigenvalue, 
where n(“VN + o) and l denote thu circumferential 
mode number and the radial mode number ,. respective- 
ly, of 0 pressure wave. 


Hie factor O n 1 ia cither a positive real or 
a positive imaginary number. Evidently the modes 
corresponding to real Op $ arc the cutoff waves 
that decay in the far- field, while the modes asso- 
ciated with imaginary values of O ni j ore the cut- 
on waves that propagate without decaying, flic 
critical valuca of v and l depend upon the 
axial Mach number M a » the rotational speed of 
rotor cap, the hub-tip ratio U and the circum- 
ferential mode number of the inflow distortion q, 
but not on the radial mode number of inflow dis- 
tortion p. 


which Is equivalent to 

DO 

%,p(P.6) - 2 C n,l ( *’ q ’ p) 'V k n,i ,p > 

1=0 


- P 2 A(p) £ p n ,z<^,P) 'MW’ 13 ) 


2-0 




Equation (9) expresses the fluctuating pressure os 
a aupcr-pouillon of on infinite number of pressure 
modes which arc characterized by the circumferen- 
tial mode number n(=vN + a) and the radial mode 
number 1. Expression (9) is convenient tor a 
modal analysis of the pressure field, but It is in- 
convenient for evaluating the near field pressure 
for two rconons. First, the series expansion with 
respect to the circumferential wave number n is 
nonuniformly convergent at the blade surface. 
Second, it involves too muny blade force coeffi- 
cients G n i<5,q,p) and Fn, i(( i4»P) (J a 
2, , . ,,v = 0, 1,2, . . . ) which are also func- 
tions of the axial position £. These difflculti- 
tlcfi, which may be crucial, especially in n case 
where the unsteady lift force distribution 
AP (p,a) cannot be specified n priori, were elim- 
inated by the method developed by Nnmba, 2 The 
method is founded upon approximating RaC^n, l‘ r ) 
by a finite series of the R 0 fc '' c ^ orm 




r fit *r) = \ ’ B . R (k , *r). 
« n,Z /_j n,l,k o o,k 

k=0 


( 12 ) 


Equation (ft) can be written in the Fouricr- 
BUssel double series from as follows; 


The number L of retained terms should be appro- 
priately large. With this method, equation (4) can 
be rearranged to obtain 


) 
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w e (r.Ov^.t) - < a -W a £ A (t)«e’ lq(0 "V > (16-1) 
8 * 


DO <0 


« C -W T' V B ' R (k t) 

“ * q£i. ptt «l»P q q -P 

Xe -iq(0-V> ( 16 - 2 ) 


'Die kernel function W’^rtd.z/l^q) 1 * governed by 
the parameters N, u>[, h, q, Wa f d and L, 

Kl’k can be resolved Lnto three parts! Kl’ 
(pronogating port), Kp(|) (singular part), and 
KP(|) (regular part). 


where q and p arc respectively the circumfer- 
ential mode number and radial mode number. t a de- 
notes a small quantity which is the ratio of the 
external fluctuating velocity to the basic flow 
velocity. Since the upwash component of external 
fluctuating velocity is given by 


Fluctuating Velocity induced by Retor Blade Row 

The fluctuating velocity q(r,tS,z) c lu)t can 
be obtained by integrating the linearized Euler's 
equation of motion as follows; 


q r,w ' ’ w q, 0 , ‘V r /V 1 + U5 T r2 (17 > 

then, one can express the upwash component of the 
extcrnel fluctuating velocity on a blade surface 
0 * WpZ by 


0(r,0,z) 


1 




c t<#B |W(r,tJ,*)]ds 

0-0 -Kyi, 


( 14 ) 


where P 0 denotea the fluid density in tile un- 
disturbed state and Q a W a ’ r^ is the velo- 

city of tho undisturbed fluid relative to the 
blades. Combining equations (ft) a nd (14) leads to 
the expression of tho fluctuating velocity in terms 
of the blade force mode coefficients, then the up- 
wosb component q(r,.>,z) e ixc can be expressed in 
the following form 


q T («v>i*> 


1 j_ 
p 0 Q ftir 


x 




L-l 

£ C o,k (f * q ' p> 


k=0 


•KT k (r,rf,z/^,q) d(. (15) 

The upwash kernel function contains parameters of 
N, h, q, uij- , 0 , and L. The detailed expression 
for the kernel function KT^ are given in Appendix 
X. As shown, KTfc is dLvlded into three parts the 
sume ns 1 T^; the propagating part KTl£), the 
singular part KT'jp , .and the regular port KTljp. 
The singular part KT'jp possesses singularities at 
, the lifting surface. 




f)=UL,Z 
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or - V »C T' 
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l + s^r 


(18) 


One easily secs that the reduced frequency uC Q 
and the interblade phase angle 2 ttS/N of Inflow 
distortion sensed by the rotor blndus are given by 

uj.Cj ts q'C^'tdj, and a a mN - q (19) 

where the integer m is chosen so that 1 2rr u/N[<tr 
is satisfied. It follows from equation ( 16 ) thot 
the Fourier coefficient Uq can be determined 
from the axial distortion velocity w e by the 
relation ’ 


, f l ». „fr.» 


e lqt} dr do 


( 20 ) 


It is convenient to suppose that the distortion 
velocity can be expressed ns the product^ 

w „(r,3) 

«F.W. 0 j( 8 ) ( 21 ) 

a n 1 


then 


where 


q.p 


a 

q 


• b 

p 


( 22 ) 


Inflow Distortion 


In this paper the inflow distortion is assumed 
to have only an axial velocity component (see Fig. 
2). If a Fourier-UesseL analysis of the arbitrnry 
shapes of inflow distortion is carried out, the 
following axial component of external fluctuating 
velocity is obtained. 


and 
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q 



e’ lq§ . 0^(5) d§ 
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r 'V k q,p' r) ‘ F l Cr) dr 


(23) 


(24) 
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When the Inflow distortion can be represented by 
N emission profiles (for N “ 1, 2, » . .) in the 
circumferential direction. 


!p, - (r - l)ir/M» 

1 (r - 1, 2, .... M) (31) 
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Wo find that 


2 2 2 

«-JjE 1 


Jyi7*N*0 2 >o -jj- q.e 

0 otherwise for q « 0, ±N, ±2N, . . 


Re 


ot£ (fir 2 + ^a) 


(26) 


Detormlnntton of Acoustic Dipole Distribution 

In the previous sections, the fluctuating 
pressure and velocity induced by the rotor blades 
have been represented In terms of the acoustic di- 
pole distribution, namely, the unsteady lifting 
pressure .'l’(riz) e^ 1 ^ on the blades. The acous- 
tic dipole distribution must bu determined so that 
the normal component ot the relative fluid velocity 
should vanish at the blade surfaces. Therefore, 
thu upwash component < l the external fluctuating 
velocity q^ w must bo cancelled by the Induced 
upwash velocity q T (iy«,a) e lu;t at the blade sur- 
faces} that is, for jZj < C a /2 , 


bet 


q^r.yys,*) e + [q lw ] 


' 0-1 


V 


(27) 


M Is the number of chordwlse acoustic dipole points 
and is also the chordwlse boundary points on a 
blade G_ u(C»qiP) are represented by 

O jN 

e 0|k ft.1.P) - C^Lq.p)/^ .In ¥j) (32) 

The quantity C^fapq.p) is calculated from equa- 
tion (23) by a collocation method, The solution is 
assumed to satisfy a Kuttu condition at the trail- 
ing edge l.c. , C* i 4 (tr*q,p) * 0. 


Pave Tone Acoustic Powe r 

The pure tono fan noise Is composed of a finite 
number of propagating pressure modes ubich depend 
on the parameters v, 2, q, .n^, M fl and h, Using 
equation (3), we can obtain the expression for 
sound pressure 


n 2 

R (k ,.r). 

n ii| l 

X ur t (n,l.O-. 1<H * )M r t (M > 

Here the suffixes + and - doooto the sound waves 
propagating forward (downstream) and backward (up- 
stream) respectively, and lU‘i(n,l,q) denotes the 
uondlmenslonnl pressure amplitude in the (n,l) 
acoustic mode, which is given in the noncoiupnce 
source case by 



Combining equation (2?) with equations (IS) and 
(18), we obtain an Integral equation tor the un- 
known acoustic dipole C 0 !<(£,, q,p) la the Conn 
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(28) 


where 


\ “ 0 


2 , , 


(29) 


We adopt the following method to cbtoin acoustic 
•'dipole distribution from equation (28). The coor- 
dinates £ nnd « are transformed into angle 
variables < p ^ and q >- such that 


i “ -(b ( /2).cos «P L 

z ■ -(C n /2)-cos q> 2 (30) 


(36-1) 


while In the case of n compact source 



Further a^., given by 


where 


5 



where 


M 2 M 2 

a ± “ “f < nW T + “> * i n n,ll “ ~i (» + 1>"*r * l n n,ll 




(35) 

denotes the axial, wave number per length r£. 

Corresponding to each circumferential mode • 
number of Inflow distortion) the circumferential 
acoustic mode numbers n (namely v) which con pro- 
pagate ore decided with equation (8). Hie pure 
tone harmonic number corresponding to each mode 
number n is not v f but v+m, where m is ob- 
tained from equation (19). The inflow distortion 
radial mode number p has no influence on whothor 
or not the fluctuating pressure caused by the in- 
flow distortion-rotor interaction can propagate, 
but it affects the amplitude of pure tone fan noise. 
Therefore, we need to sum all sound pressures cor- 
responding to each circumferential mode number q 
and radial mode number p of the inflow distortion, 
in order to obtain the pure tone fan noise caused 
by the arbitrary shaped inflow distortion. 


- |KP ;t <n,Z.q)| 3 2 (nu> T + a))|n n)l j i f(a i + ntq r + to) 2 (40) 

Et4(n,l,q) is the modal component of the dimension- 
less acoustic power. 


Two- Dimension "l Theory (Strip Theory) and Quasi 
Threc-Dlmensl in l Theory 

A two-dimensional cascade calculation Is ob- 
tained as the limit of large cylinder radii. In the 
notation of this paper, this limit corresponds to 
h - 1 and N — « simultaneously, while keeping 

27r/(NC a Vi+«4) » t constant at each radiuB of the 
annular blade row and q/N = constant, where t is 
the pitch-chord ratio and q/N is the'ratlo of in- 
flow distortion circumferential wave length to 
cascade pitch. 


From equation (15), we can get the expression 
of the propagating axial fluctuating velocity in 
the following form, 


w 


a,p 


,(r,0,z) “ 


X c t ( n0+a ± z ) .HP ± (n, l,q) C l <v-*«)1«V (36) 



Following Morfey's approach, Z8 the axial component 
of tile sound intensity lJ± of pure tone fan noise 
(harmonic number | = l, 2, 3, . , .) referred to 
the duct-fixed coordinate syutem is obtained in the 
form 




(i + H 2 ) <P J 'vr>) +(w 2 /p o a 2 ) P 0 W fl (wi 2 ) 


(37) 



c ■ptf'/R 2 (k 
a o o n' n. 


I'O (nuij. + u>) p 2 


< . 

i X |0 nj J|np ± (n > l,q)| 2 / /8[a ± + (n(n T + m)] 2 (38) 

where ( ) denotes iihs time average value. Also, 
of course, p and W in equation stand for the 
''real parts of the original complex forms of p and 
W.' 


' ,, Tlicrefore, the dimensionless acoustic power 
E* with respect to each pure tone fan harmonic 
(for example j = +m=±l corresponds to 1st harmon- 
ic) is obtained by 




P wVr * 2 ' 

'o a a T / 


“2 

EE 

q n=n^ 


V EI^(n,l,q) 
1=0 


( 39 ) 


In the quasi three-dimensional method, first, 
the chordwlse unsteady force distributions (nsmely, 
acoustic dipole chordwlse distribution) are calcu- 
lated at L numbers of radial positions from hub 
to tip using the two dimensional cascade calcula- 
tion. Second, those dipole distributions on a 
blade are resolved into Fourier-Ccssel form with 
equation (11) and then the pure tone fan noise and 
modal structure in three-dimensional annular cas- 
cade are computed from equations (39) and (40). 


Numerical Results 


The numerical calculations are carried out In 
this paper for two different rotor configurations. 
Their geometric parameters are given in Table I. 

The fan denoted case No. I has a high tip-speed and 
a large number of bLadcs, while the No. 2 case fen 
has low tip-speed and a smalL blade number. Under 
those conditions, the pitch-chord ratio S, the 
stagger angle f, and the ratio of the relative 
flow-velocity to the axial flow velocity Q/W a 
vary along the span as shown In Fig. 3. Reduced 
frequerfeics of blade force on the blade due to an 
inflow distortion in cases No. I and No. 2 are cal- 
culated an Kj^ = 0.136 q and k 2 = 0.272 q, re- 
spectively. The pure tone acoustic power has been 
computed for different combinations of the circum- 
ferential mode number q and the radial distribu- 
tion of the inflow distortion. In order to clarify 
the three-dimensional character and the accuracy of 
available two-dimensional theory, the results ore 
compared with those of the strip theory method and 
the quasi three-dimensional method. 

The finite series approximation of L = 7 and 
M = 6 were used for the calculations in this paper. 
In accordance with the increase of circumferential 
wave number q, the number of propagating modes 
(n,l) increases, and therefore the accuracy of the 
calculation of the KT(|j) term in equation (15) 
seems to decrease. In the three-dimensional cal- 
culation, the degree of accuracy seems to be great- 
er than in the two-dimensional calculation, because 
the three-dimensional calculation has three- 
dimensional effects in the radial direction. So 
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the calculation fat* a large number of propagating 
acoustic modes will require an Increase in M. 

Acoustic Power Variation with Circumferential Mode 
Humber of Inflow Distortion 

In Fig, 4, the dimensionless total acoustic 
powers from the 3-D calculation of the upstream and 
downstream propagating sound waves are plotted 
against the circumferential mode number q of the 
Inflow distortion. The coefficient Aq in equa- 
tion (16-1) was hold constant over all values of 
q .tnd the radius. In this section, we diucuss the 
n-p.’i rical resuLts for the noncompact source case. 

The din -tslonicss acoustic power at q » 10 agrees 
with tl'iw Namba'a results. 1 However, acoustic power 
levels or, other values of q are different from 
Nambo's results, because in Numba's calculation thl' 
i reduced frequency Is assumed ns constant ultheuc 
jregard to the variation of q, In the upper part 
of Fig. 4, the variation oi the generated acoustic 
modes with the inflow distortion mode number q is 
shown, in the notation (vj,’) Vi corresponds to 
the harmonic number of the pure tone and I corre- 
sponds to the radial mode number of the acoustic 
•mode. Therefore the circumferential mode number of 
the acoustic mode n is calculated with the equa- 
tion n ■ viN-q whore N is the blade number and 
q 1 b the circumferential mode number of the Inflow 
distortion. 

I 

; Figure 4 shows that the number of acoustic 
modes Increases with an increase in q. The acous- 
tic powers of both upstream and downstream propa- 
gating sound waves are largest at about q “ 30, 

This region of large acoustic power corresponds to 
n large number of 1st harmonic acoustic nvidea. in 
Fig. 4, it it< sueu that the total acoustic power of 
downstream propagating waves is higher than that 
upaivennt propagating waves in tile whole range of q. 

Figure 5 shows the variation of fundamental 
pure tope component with circumferential wave num- 
ber q of ini low distortion for case Ho. 1. For 
values of q less than 20, the fundamental tone 
acoustic power level of downstream propagating waves 
is higher than that ot upstream propagating waves, 
but for values ol greater than 20, this tendency is 
reversed. Total acoustic power and the acoustic 
power of Hu* fundamental pure tone Cor case No. 2 
ore shown In Figs.. 0 mid 7, The total ncottscic 
power radLnted downstream la higher than that radi- 
ated Upstream, while for the fundamental pure tone 
the downstream acoustic power is somutlsiea lower 
mid sometimes higher than the upstream acoustic 
power. 


Compact and Noncompnct Sources 

j 

1-k In Figs. 4 to 7 , comparisons arc made between 
| the compact source prediction based on distributed 
sources only along a line in the spnnwi.se direction 
‘(compnct source) and a rigorous prediction based on 
distributed sources along both the airfoil chord 
‘and span (nencoropnct) for the upstream and down- 
stream propagating sound waves. The divergence of 
the compnct and noncempact source predictions in 
both the upstream and downstream propagation cases 
is considerable, especially in the lawer reduced 
■ frequency range for small values of q. This ten- 
'dcncy is similar to that observed by Kajl 12 in a 
two-dimensional cascade. The maximum acoustic 

- 


power difference is 15 dn In tho upstream case and 
6 dB In the downstream e-n-fv, Tlio degree of differ- 
ence decreases as the reduced frequency of tho 
blade forces Increases, 

At the higher reduced frequencies , tho unsteady 
forces on the blade are concentrated near the load- 
ing edge and so the phase variation of source in 
the chordwluo direction is small, In addition, in 
the downstream direction, tho wnva length of sound 
Is stretched due to fluid flow and therefore the 
ratio of wave length of sound to the chord length, 
a retarded time effect, Is smaller than in the up- 
stream ease. Therefore, the difference in acoustic 
power between noncempact and compact suurcos is less 
noticeable In downstream case (sec Fig. 5). 

In tho case of the fundamental tone (Fig. 5), 
the downstream ucoustic power predicted by the com- 
pact source model la higher by n maximum of 7 dll 
than that of the noncompnct source modeL and they 
approach one another in rgo higher source reduced 
frequency range. In the upB cream cuse^ the ncous- 
tic power level of the eompoet source prediction is 
considerably higher than that of the noncempact 
source prediction in the low reduced frequency 
range, but in the higher reduced frequency range, 
the ocuuscic power of the latter Ls higher than 
that uf the former by several dll. 

Figures 6 and 7 for case No. 2 also show that 
the acoustic power difference between compact and 
noncompact sources is smaller than In case No. 1, 
but still appreciable. 


Acoustic Power Variation Duo to Inflow Distortion 
Radial Distribution 

Figures 8(d) and (b) sin ’he acoustic power 
variation ol the fuiidament.il for upstream and 

downstream propagation. Tin , relations were for 

ease No. 1 using Iho 3-1) modi ;. In these figures, 
acoustic powers fur ') cases ol inflow distortion 
radial distribution are shown. The 5 canes arc lor 
Inflow distortions with radial extents of 8, 25, 

42, i»7 and 100 percent uf the a par from starting 
from tlio blade tip (see sketch in fig. 8), these 
figures show that at tow values oi q, nil oi the 
upstream acoustic puwor is generated by t\U* outer 
25 percent of blade- In flow distortion Interaction, 
At higher values ot q, the 25 percent of blade- 
inflow distortion Interaction still accounts for 
over hall of the acoustic power. For the down- 
stream wove a ulmLlar behavior was observed that 
involved tho 42 percent of blade-inflow distortion 
interaction. 


Hodnl Structure 

Figures 9 and 10 show the variation of radial 
mode level for the fundamental pure tone noise for 
various values of q for both enses No. 1 and No. 

2. It Is scon that tho aoiuutic power in the first 
radial mode Is generally Higher Chan that in Che 
second radio L mode. 

In the case of upstream propagation, for case 
No. 1 (Fig. 9(a)) the acoustic power of the first 
radial roodu is higher than second rodinl mode, by 
more than 7 dll for q less than 40, In case No. 2, 
(Fig, 10(a)) the second radiul mode level is close 


to the (trot radial mode at q • U and 16, however, 
for •>! othar valuaa of q, tha fonaar la lower by 
about S dB than tha lattar. 

In tha caaa of downstream propagation (Tig. 
9(b)) caaa No. I shows that tha aacond radial mode 
l aval la lowar by nor a than 3 dB than first radial 
aoda laval. Othar radial more lavals ara lowar yat 
than aacond radial mod* laval. In caaa No. 2 
(Pig. 10(b)) tha aacond radial mod a loval la lowar 
than tha flrat radial component laval aacapt In tha 
ranga of high raducad frequency. Tharafora, In 
almost tha antlra range of q and for both up^ 
•traaa and downstraam caaas, It a cams that only the 
first radial mode might ba enough for acoustic 
liner designs If auppreaaton requirements ara mod- 
ast. 


Accuracy of Two-Dimensional Method for Pure Tone 
Pan Noise 1'rcdictton 

Comparisons ara made of tha acoustic powers 
calculated from the three-dimensional method (3-D), 
the two-dimensional method (2-D) and the quasi 
three-dimensional method (quasi 3-D), In order to 
clarify the three-dimensional effects and the accu- 
racy of the two-dimensional method for fan noise 
prediction. Figure 11 is for case No. I and Fig. 

12 la for case No. 2. 

Numerical results for both upstream and down- 
stream total acoustic power for case No. 1 (Figs. 
11(a) and (b)) show that the acoustic power differ- 
ence between the two-dimensional and three- 
dimensional calculation is less than 3 dB. The 
quasi three dimensional calculation results based 
upon the two-dimensional unsteady force are as 
close to the 3-D results as are the 2-D results. 

Reaulta for the fundamental pure lone (Figs. 
11(c) and (d)) show that the level differences be- 
tween 2-D (also quasi 3-D) and 3-D ate smaller than 
for the total acoustic power, litis is true for 
both the upstream and downstream waves. 

For case No. 2, Fig. 12(a) shows the upstream 
total acoustic power difference between >-D and 2-D 
to be about 3 dB at some values of q, but the 
acoustic power difference between 3-D and quaai 3-D 
is less than 4 dB. For the downstream wave (Fig. 
12(b)), the total acoustic power difference is 
generally 1 drt or less in the region between q * $ 
and q *• 24. Figures 12(c) and (d) show similar 
close agreements for the fundamental tone. 

The survey based on the above data Indicates 
that the quasi 3-D method may be used for pure tone 
fan noise prediction with an acoustic power esti- 
mation error generally lass than »2 dB. 


Comparison Medal Structure Between 3-D and Quasi 
3-D Mathilda 

Figure 13 niiuws ct\e comparison between modal 
structure of the fundamental pure tone baaed on 3-D 
and that based on quasi J-D methods. Figures 13(a) 
and (b) are for case No. 1 and show the upstream 
and downetream cases, respectively, while, Figs. 
13(c) and (d) are for case No. 2. In both cases. 
No. 1 and No. 2, the agreement of acoustic powers 
of corresponding acoustic modes between 3-D and 
quaai 3-D methods la good. In both rotor cases, 


tha acoustic power dlffaranca of corresponding 
acoustic modes la lass than 3 dB, and In tha region 
of low reduced frequency, the power difference 
approeches aero. In accordance with these numerical 
reaulta, tha quasi 3-D method also seems to be ade- 
quate (or acoustic modal structure prediction In 
the three-dlmen* tonal annular fan •• Inflow distor- 
tion interaction problem. 


Joncluding Remarks 

Comparisons ara made among tha predictions 
basad on thraa-dlmanslonal theory, atrip thaory and 
quasi thrae-dlmanalonal thaory. In order to clarify 
the accuracy of aval labia two-dlmenalonal thaory 
for the prediction of pure tone fen noise due to 
the Interaction of Inflow dtatortlon with a subsonic 
annular blade row. The theories wera derived with 
the eld of an unsteady thraa-dlmanslonal lifting 
aurfaca theory developed by Nambe.l Both compact 
and noncompact sources were considered. Several 
numerical calculations were carried out and tha 
following results were obtained. 

1. The compact source prediction in a three- 
dimensional annular blade row cen overestimate the 
upstream radiation by 15 dB and the downstream radi- 
ation by 6 dB. The degree of overeatlmatlcn depends 
on blade row characteristics and reduced frequency 
end interblade phase due to inflow distortion. 

This result implies that noncompact source affects 
must be considered. 

2. The upstream radiation seems to depend 
heavily on the tip region ot bl.’ je- ini low distortion 
interaction for both upstream and downstream cases. 
The extent of apanwise contribution to the sound 
was greater in the downstream case than in the up- 
stream case. 

3. The numerical results show that for acoustic 
liner design the consideration of only the first 
radial acoustic mode is enough when only mod os t 
suppression is needed. 

4. Strip theory and quasi three-dimensional 
theory ’»ased on the two-dimensional calculation of 
unsteady forces are reasonably adequate for tan 
noise prediction. The quasi three-dimensional 
method Is more accurate for acoustic power and for 
modal structure prediction, with an acoustic power 
estimation error of *2 dB or less. 
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Appendix - Kernel Function -Hk 

The kernel function KT^ la composed of the 
three parts KT^. KT'f), and KT'f', which are 
given es follows: 
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TABLE I. - FAM PARAMETERS 


i 

Case 1 

Case 2 

Rotor blade number, N 
Axial chord length. 

40 

15 

dimensionless , C~ 
Rotor spccd/oxtnl flow 

0.05498 

0.25836* 

speed, u>f 

Rotor relative Mach 

2,474 

0.627 

number, Mj. 

Axial flow Mach number, 

0.934 

0.865 

M- 

0.35 

0.59G 

Hub/tip ratio, h 

0.4 

0.46 


0\w 

0 ? 
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Figure 4. - Variation of total acoustic power with inflow dis- 
tortion circumferential mode number, case no. 1. (Uni- 
form radial distribution of inflow distortion. > 
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Figure 6. - Variation of total acoustic power with in- 
flow distortion circumferential mode number q, 
case no. 2. lUniform radial distribution of inflow 
distortion. I 
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Figure 7. - Variation of fundamental 
pure tone noise with inflow dis- 
tortion circumferential mode num- 
ber q, case no 2 (uniform radial 
distribution of inflow distortion!. 
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Figure 9. - Variation ol radial mode levels ol fundamental 
pure tone noise with inflow distortion circumferential 
mode number. lUniform radial distribution of inflow 
distortion. ) 
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Fiyiire 10. - Variation of radial mode 
levels of fundamental pure tone 
noise with inflow distortion cir- 
cumferential wave number. (Uni- 
form radial distribution of inflow 
distortion. I 
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Figure 11. - Comparison of acoustic powers among 3-D 
method. 2-D method, and quasi 3-D method - case 1. 
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Figure 11. - Concluded. 
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Figure l£ - Comparison of acoustic powers O mong 
3-D method. ?-D method, and quasi 3-D method - 
case Z 
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Figure 1?. - Concluded. 


30. 



la* UPSTREAM CASE. NO 1. 



circumferential mode number of infiow 
DISTORTION q 

Ibl DOWNSTREAM CASE. NO !. 

Figure 13. - Comparison of modal structure of fundamental 
pure tone noise between 3-D method and quasi 3-D method. 
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Figure 13. - Concluded. 



